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Abstract. Discussion is given of non-linear soliton behavior including coupling 
between electrostatic and electromagnetic potentials for non-relativistic, weakly 
relativistic, and fully relativistic plasmas. For plasma distribution functions 
that are independent of the canonical momenta perpendicular to the soliton 
spatial structure direction there are, in fact, no soliton behaviors allowed because 
transverse currents are zero. Dependence on the transverse canonical momenta 
is necessary. When such is the case, it is shown that the presence or absence 
of a soliton is intimately connected to the functional form assumed for the 
particle distribution functions. Except for simple situations, the coupled non- 
linear equations for the electrostatic and electromagnetic potentials would seem 
to require numerical solution procedures. Examples are given to illustrate all of 
these points for non-relativistic, weakly relativistic, and fully relativistic plasmas. 
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1. Introduction 



Since its discovery in 1959, the Weibel instability ( Weibell . Il959t iFriedl . Il959l ) has 
been, and continues to be, a subject of intense research. Part of this focus is due to 
the fact that the Weibel instability excites aperiodic fluctuations, i.e., purely growing 
waves that do not propagate out of the system. Therefore, the Weibel instability 
has been considered to be responsible for the creati on of seed rn agnetic fields in the 
early universe (Schlickeiser_and Shukla, 2003; Sakai et al. , 12004 LSchhckciscr, 2005), 
which act as a progen itor to the large- scale magnetic field that we observe today 



in all (spiral) galaxies ( Beck et al. . 19961) . Likewise, Weibel fluctuations can provide 



the necessary dissipation of bulk velocity that leads to the formation of shock waves 
(|Tautz and Lerchel l2006th . n highly relativistic jets that occur during events such as 
GRBs (Gamma-Ray Bursts) or in extreme obj e cts su ch as AGNs (Active Galactic 



Nuclei), plasma instabilities (jSchlickeiser et al.l 120021) and particularly the Weibel 



instability are, thus, ultimately responsible for the radiation observed at Earth 



Weibel modes and their asssociated non-linear structures ( Bvchenkov et al. . 20031 ) 



also play a rol e in the radiative cooling of rela tivistic particles in blazar and gamma-ray 
burst sources ( Schlickeiser and Lerche . 20071). Furthermore, Weibe l modes can also be 
excited in quantum plasmas ( Haasl . 20081 : iHaas and Lazar , 20081) . thus generalizing 
the classical Weibel instability equation; here, a connection has been made to (dark) 
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soliton waves (jShukla and Eliassonl . 12006 ). It is also wo rth noticing that there exist 



exper imental verifications of Weibel (|Wei et all 12001/021) and soliton (jBorghesi et al 



l2002h modes in laser plasmas, thus emphasizing the broad applicability of the 
underlying mechanism, which converts the free energy from anisotropic distributions 
into magnetic field energy. In analytical studies, such so l iton modes have been u sed 
to create magnetic turbulence (iKingsep and Sudan . 1973 : Weibel and Chen . 1977? ). 

The underlying analytical investigation of the non-linear aspects of the classic 
Weibel instability made use of the fact tha t the classic Weibel instab i lity e xcites 
only transverse, electromagnetic fluctuations (jSchaefer-Rolffs and Lerchd . 120061 ). For 
asymmetric distributions, it was shown that the range of unstable wavenumbers is 
reduced to one single unstable wave mode, which reminds one of solitary structures 
that are based on single wavenumbers, too. For the case of transverse electromagnetic 
modes, it was shown that, depending on the exact form of the distribution function, 
spatially limited structures are produced (solitons). 

From the radia tion pattern of particles scattered in soliton modes 
2009h . it is known that there are many similarities to syn- 



( Schaefer-Rolffs et al 



chrotron radiation. Motivated by the fact that the energy output of the particles is 
mostly referred to as synchrotron radiation and inverse Compton scattering, it is a 
necessity to take into account other radiation processes, also because synchrotron radi- 
ation requires the presence of a background magnetic field. Unt il now, the origin and 
theni e chanism fo r such a field have been disc ussed intensely ( Medvedev and Loebl . 
19991 : iMedvedevl . 120071 : iMedvedev et all . l2007t ). The train of thoughts is the fol- 



lowing: (i) in reality, particle velocity distributions should almost always be some- 
what asymmetric; (ii) in this cas e, unstable p l asma, modes are monochromatic, as has 
been demonstrated analytically ( Tautz et al.l . l2006h : (iii) such isolated structures can 
lead to soliton modes as has bee n shown for purely electromagnetic Weibel modes 
( Schaefer-Rolffs and Lerchelliooeh : (iv) in relativistic outflows such as for GRBs and 
AGN jets, all kinds of plasma instabilities are expected to arise. Hence, it is most 
important to discuss the radiation pattern for such scenarios that is generated by 
particles being scattered in such magnetic structures. In principle, the method of ob- 
taining the differential frequ ency spectrum is similar to that for synchrotron radiation 
( Schaefer-Rolffs et al. . 20091 ) . The basic difference is that for synchrotron radiation 
the magnetic field is considered to be constant and the electron moves in circles per- 
pendicular to the magnetic field, while in the case of the soliton the electrons move 
mostly linearly and are deflected via the Lorentz force. Thus, the radiation is not 
produced by acceleration through a constant background field but, instead, is caused 
by an interaction of the electrons with highly varying magnetic and electric fields. 

Inspired by the procedure shown in ISchaefer-Rolffs and Lerche (|2006[ ). a number 
of subsequent, more detailed, investigations revealed the exotic properties of the 
(linear) Weibel instabil ity that are unfolded in the case o f totall y asymmetric 
distri b ution functions (Tautz et al.l , l2006t iTautz and Lerchei 12006 d : [ Tautz et al 



20071 : iTautz and Lerchd . l2007a ). For such distributions, the electrostatic and 



electromagnetic wave modes are coupled, and it was shown that any unstable Weibel 
mode must be isolated, i.e., restricted to a single discrete wavenumber. Specific 
examples for the distribution function illustrated that isolated Weibel modes are 
excited. Even i f one allows for a small re al part of the frequency, the isolated Weibel 
modes persist ( Tautz and Lerchd . 2007i). Such weakly propagating unstable modes 
are important for oblique wav e propagation, because for such cases the growth rate 
of unstable waves is maximal ( Bret et al. ■ 12004 120051 ). It is then appropriate to ask 
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how non-linear soliton modes are influenced when one includes the coupling between 
electrostatic and electromagnetic potentials. The purpose of this paper is to explore 
that question. 

The paper is organized as follows. In Sec. [21 the Vlasov equation is transformed 
and the non-linear wave equations are derived. In Sec. |3l the non-relativistic soliton 
behavior is derived and examples are given for two basic types of distribution functions. 
In Sec. m the weakly and the fully relativistic behaviors are derived, which requires 
approximations as regards the transformation of the volume element in momentum 
space. Furthermore, in Sees. [3] and SI solutions are derived for two limiting cases 
regarding the relative values of the electrostatic and vector potentials. In Sec. [Sj 
the results are summarized and discussed. Detailed explanations of the integral 
transformations in the cases of non-relativistic and relativistic plasmas are given in 
[Appendix A| and [Appendix B[ respectively. 



2. The Relativistic Vlasov Equation 



Throughout the derivation of the non-linear wave equation and the instability 
conditions, the vector potential, A ~ {A^, Ay, A^), and the scalar potential, $, will 
be used, with 



IdA 
B^V xA. 



- V<I> (la) 

(16) 

For a soliton wave, observed to be moving at a constant velocity /3c, the easiest way to 
handle the non-linear equations is to transform to a reference frame moving with the 
soliton so that, in such a reference frame, the potentials A and <& are functions solely 
of the spatial coordinate, x, and are independent of time. The electric and magnetic 
fields are then given by 



E 
B 



(0, ~A'^, A' ), 



{2a) 
(26) 



where primes denote differentiation with respect to x. The Vlasov equation can then 
be expressed a^ 



df e 
ox c 



-VyA'y 



(3) 



a/ 

dpx '^Fy 

The characteristics of Eq. with respec t to the y and z coordinates are giv en through 
dp()^,z)/dj4(j, 2) = — e/c (compare with ISchaefer-R olfFs and Lerchel 2006 ). For the 
solution of the characteristic equations introduce i37(y 2) through 



^y,z) - P(y,z) + 



'-A 

C 



(4) 



Then, Eq. ([3]) simplifies considerably, because '^(y^z) ^-re constants and thus the partial 
differentiations of the distribution function / with respect to Wy and Wz vanish, 
yielding 

df e dj 



= 



rwy dx c dpx 



1 / 

Wy 



A. A' 



-(' 

7717 v 



Az A' 



(5) 



X While our notation differs from that used in classical relativistic mechanics, it is traditional to the 
field of plasma physics. 
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with / = f{x,Px, T^y, TUz) and where 



7 



1 



1 



c 



(6) 



It is instructive to first consider the non-relativistic hmit [set 7 = 1 in Eq. ([5])] 
because this hmit aUows exact solution of the Vlasov equation for arbitrary $, Aj,, 
and Az- The non-linear coupling of the electrostatic and electromagnetic potentials 
in the Maxwell equations can then be investigated simply in the weak coupling limit, 
where "weak" will be defined later. 

Thereafter one can consider the relativistic particle behavior where, as will be 
shown, coupling of the electrostatic and electromagnetic potentials is considerably 
more involved, although the basic procedure from the non-relativistic situation can be 
used, mutatis mutandis. 

3. The Non-Relativistic Limit 



In the case where one deals solely with non-relativistic particles one can set 7 = 1 in 
Eq. (O. Then one has 



TO dx dpx 



+ — fujy-- Ay) A' + — (mz--Az] A'^ 
m,c \ c / " TOC V c 



Set p^ — u m Eq. Q to obtain 

^ df „ df 
= + 2eTO 
ox ou 



mc V 



Wy Ay]A 



" m,c \ 



With 



w = u + 2em$ 



ZUy — Ay 



wz- - Az) Ai 
m,c \ c 



-A, 

c 



(7) 

(8) 
(9) 



Eq. ([8]) provides the general exact solution / = f{w,vjy,Wz) by direct substitution|§| 
Consider then Maxwell's equations in the limit of non-relativistic particles. For 
a plasma with different constituents — denoted by the index a with, e.g., a = e for 
electrons and a = p for protons — one has 



4>" = -47r^e„ /d^p/, 

-47rV-^ / d?PP(y.z)fa 
„ TOaC J 



A" 



(10a) 
(106) 



In terms of Wy^ and Wz the Maxwell equations take the form (see [Appendix A| ) 



{•OO pOO pOO Q 

^''^SttV^Cq / d^ / AuJy / dWz ^/l^ !a{Wi, + S,,Wy,Wz) 
Jo J -00 J -QO C*? 



(11a) 



and 



pOO pOO pOO 

AI^SttY,— dwy dwz 

TOq,c Jo j -00 J -00 



(116) 



§ Consider Eq. ^ in the form df /dx — {dg/dx)(df /du) = 0. Following the theory of linear partial 
differential equations, set g{x) as a new variable. Then df/dg — df/du = 0, which has the general 
solution / = f{g + u). With w = g + u, this is precisely as given. 
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with A_L — (0, Ay,Az), 'c^± = (0, zuy, Wz) and w = ■w^, + ^, where 



— A_L 

c 



(12) 



Solution of the coupled electrostatic potential, Eq. ()llap , and the electromagnetic 
potential, Eq. (jll6p , depends sensitively on the assignment prescribed for the particle 
distribution function f{w,Wy,Wz). Note that if f{w,Wy,7jJz) is a function solely of 
w and not of zUy,'CDz explicitly, then A± — 0. Then only $ is eligible for a soliton 
structure. To illustrate this point more closely, consider two cases, where f(w, Wy, zuz) 
takes on the functional forms 



or 



/ = fa.O CXp ( — 



fa = fa,o exp ( ) exp I 



(13a) 



(136) 



with fa.O, Wa, and Wa constants. 

Direct integration of the right-hand sides at Eqs. (|llap and (jll6[) is possible with 
the expressions from Eqs. (|13ap and (|135p . The results are in the case of Eq. (|13ap 

2eam„ ^\ ^^^^^ 

(14&) 



$" = -47r5/2^/,^oeai«;^/'exp 



A'[ =0 



and in the case of Eq. (|13&l 



3/2 2 



exp 



2 C 7T^ 



Wo 



together with 



4 = 4-^/^E/'^^o(t) 



A? 



Wo 



(15a) 



■A, 



X exp 



2 ^2 + 

y -2 



(15&) 



with i € {VtZ}. Note that as Wa ^ oo Eqs. (|15a|) and (jl5&P reduce to Eqs. (I14ap and 
(|146|) . respectively, as required. The determination of a soliton structure in both cases 
can then be readily given for different plasma conditions. 

For instance an electron-positron plasma with identical plasma characteristics 
(i. e., fafl the same for both species) allows one to write Eq. (|14al) in the form 



„3/2 



which integrates once directly to yield 



exp 



2 6qi TTliQ^ 



(<i>')^ = i^^/.coshf^ 



exp 



const. 



With vj/ = 2mQeQ$/ii;c( one has 

,2 47r5/2u;5/2 /4TO2g2 



(*') 



Too 



/q cosh^' + const. 



(16) 



(17) 



(18) 
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Figure 1. The behavior of Ay as a function of the normahzed variable, f , for 
positive and negative values of Bq at A = j4«, respectively. 



Set the constant to a negative value, i.e., — cosh^'* (a positive constant would 
automatically keep ^' growing and so cannot represent a soliton), yielding 

(*')^ = IGn^^^wl/^famael (cosh^' - cosh**) . (19) 

Then if 4' exceeds it does so thereafter for all coordinates and so does not represent 
a soliton. If ^E* is less than '^|^, then Eq. p9)) is not valid because it would yield (^I^')^ 
negative. (For negative ^E*, the argument can be reversed.) Hence the only solution is 
a spatially unbounded potential ^E* which does not correspond to a soliton. 

Such is in the line with the small amplitude limit that provides a dispersion 
relation of the form ~ ujp + k^v^-^ which, for uj — 0, indicates k — zLiujp/ vth 
and so represents spatially unbounded growing or decaying modes. Such unbounded 
structures do not represent soliton modes and so are to be discarded. The only solution 
left for Eq. ^ is then $ = 0. 

In the case of Eq. (jl36|) there is a common factor in Eq. (|15a|) of 

exp - I — . 2 

SO that for the electron-positron plasma one again has <E>" > everywhere and so 
no electrostatic soliton is available. Here, "common" means that the factor does not 
depend on summation over species and so can be brought outside the summation. 
Thus $ = is the only acceptable solution of Eq. (|15ap . Equation (jl5&[) takes on the 
generic form 

A'^ = P Aj exp[-a {Al + Al)] (20) 
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t 


\^y\ 




c 




> 



Figure 2. The behavior of the absolute value oi A'y as a function of the normalized 
variable, Q. 



with a, /3 > and j S {y, z}. Then either A'' /Ay = A'/JA^ or one of the components 



Ay or A- 



is zero. Consider that A, — 0. Then 



which integrates once to provide 



— [const ■ 
a '■ 



Let Ay have an extremum on Ay 



exp{~aAl)] . 

Ai,. Then write Eq. in the form 



(4)' = ^o{l-expha(A^~A^)]} 



(21) 



(22) 



(23) 



with 



-aAi 



exp^^ 

so that A'y = Q oil Ay = Ai, . Eq. ([25]) shows that has the structure given in Fig. 1 
with the two solutions only just touching on = but not crossing because 

Al^QonAy = A,, for A, ^ 0, 

representing both a background magnetic field Bq together with a soliton pulse 

is 



superposed. The absolute value, \A'y\, is illustrated in Fig. O Although Eq. 
not analytically resolvable, the structure of the solution is indeed seen by considering 
turning points and roots. 

The point, then, is that exact non-linear solutions to the coupled Maxwell- Vlasov 
equations can be arranged to provide a rich variety of structural behaviors depending 
on the specific choices made for the distribution functions for the particles. In order 
that there be an electrostatic soliton component in Eq. (|llap and (jllSp one requires 
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that the precise symmetry invoked between electrons and positrons be broken. Such 
a break can be arranged in one of two ways: either one chooses different distribution 
functions for the positive and negative charged components or one chooses different 
particle masses for the positive ions with respect to the negatively charged particles. 

While it was shown that one class of solution is that when one of the components 
Aj (with j g {y,z}) is zero, that does not rule out other classes of solution where 
Aj is non-zero. In addition one cannot provide a general statemenet on the field 
direction relative to the background field because each choice of distribution function 
must be investigated for its effects on the non-linear system. It also by no means 
follows that a plasma with no background field cannot posses solitons of electrostatic 
or electromagnetic or coupled properties — again, each distribution function chosen for 
each species provides its own properties- which was part of the aim of the special 
illustrations chosen. 

For example, if one were to treat with finite electrons and infinitely massive ions 
then, with the distribution function from Eq. (|13a|) and ()13bl) . one obtains 



3/2 2 

= 47r5/2/e,o ; 1 exp 

We + Wi 



2em ^e\^Al + Ai 



Anen, (24a) 



where n is the ion number density and 



A'j=4n^/^fe^,{^ \ "^' 4, exp 



2em ^ /e\2 Al + A 



(24&) 



with j e {y,z}. 

Note that one class of solutions is when Aj — and when such is the case Eq. (|24aP 
has the structural form 

= b exp(*) -gn, 6 > 0, .g > (25) 

Eq. (|25|) integrates once immediately to yield 

(*')^ 26e* - 2571* + A, (26) 

where A is a constant. Consider now in detail the structure of Eq. ([26]) . Write 

= ^' — A/ {2gn) where one has 

(i^' f ^ Be" - 2gnv (27) 

with B = 2feexp[A/ {2gn)]. Further simplification occurs when one sets x = {2gn)^^^^( 
and B/{2gn) — when one has 



^^j =R'o''-^. (28) 

Note that A, the constant of integration, appears through the relationship between B 
and R. Only for specific ranges of A can one expect to obtain a soliton behavior as 
we now show. On v — Q, Eq. (|28|) yields 

dv 

SO that two branching structures exist. If a soliton is to exist then it is necessary that 
dv/d^ — somewhere, implying a value i^* such that 

R^e"*^!^^ (29) 
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Solutions to Eq. (I^Hl) exist only when < and, under that condition, there 
are two positive values of v-i,, namely, vi^ and vu, where i>\j > vi^ without loss of 
generality. Now Eq. ([^5)) requires (di^/dC)^ ^ everywhere. One therefore has to 
distinguish three cases: 

• l\\vi^<v< v\j one has R^e'^ < v so that Eq. cannot be satisfied. Hence the 
potential regimes for solitons are either v > i/jj or v < i>i^. 

• For f > lAj one has (dz//dC)'^ > so that no bounded structure is possible. 

• Equally for v < vj^ one has the same argument with the added problem that 
(di//dC)^ grows indefinitely in v < 0. Thus no electrostatic soliton is possible. 

The point of this illustration is that the non-linear equation describing potential 
solitons is completely changed from that for an electron-positron plasma. The electron- 
positron plasma illustrated earlier has no turning points for the electrostatic potential 
taken on its own and so there are no soliton structures. In the case of the ion-electron 
plasma the corresponding field equation admits of two turning points and the analysis 
of the field equation in each domain has to be undertaken. It is this fact that the 
example has been used to illustrate. Every choice of distribution function must be 
investigated anew for the turning points and the corresponding analysis, i. e., the above 
case-by-case investigation, has to be undertaken in each regime. 

If one chooses not to have one of the components Aj = (with j G {y, z}) then the 
electrostatic and electromagnetic components mix, representing a hybrid soliton. The 
point to make is that it is the symmetry breaking (for different distribution functions 
or the different masses of the charged particle species) that permits different types of 
soliton patterns, and that has been the purpose of the examples given here. 

4. The Relativistic Behavior 

In the situation where one cannot set 7 to unity, a considerably more complex behavior 
arises but, at the same time, one has the transverse {vy,Vz) particle velocities limited 
to ±c unlike in the non-relativistic situation where the corresponding limits are set as 
±00. This fundamental difference in behavior alters radically the soliton character. 
Here one has 




(30) 



where 




so that 




(32) 



The characteristic is given through 




(33) 



which has the basic structure 




a'{x) ~-b'{x). 



(34) 
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The relativistic form of the characteristic equation is not analytically tractable in 
general. However, in the case of interest where the influence of an electrostatic 
(electromagnetic) potential on an electromagnetic (electrostatic) soliton is sought one 
has either \a'\ <C /7 or \a'\ ^ \b' \ /j. In these two situations it is possible to derive 
approximate characteristics. Also, in the weakly relativistic limit where one can write 
7 = 1 + e with only first order in e terms being held in Eq. (|33|) . it is possible to 
perform a complete analytic investigation, thereby illuminating the transition between 
non-relativistic and relativistic limits. Consider this case first. 



4-1- Weakly Relativistic Behavior 



Consider Eq. ([M]) with y = 1 + e and e considered small, i. e., |e| <^ 1. Then to first 
order in e one has 



= -a\x)-b'{x)+b'{x)e. 

ax 



Then 



da: 



;e-M 



.e-''(-)a'(x) + -f (e-"), 
da; 



with the solution 



eoe 



b(x) 



+ l-eb{^) / da;' a' {x')e-^^'''\ 



(35) 
(36) 
(37) 



where eo is a constant: the characteristic constant. Here, for a{x) one has 



a'{x) 



a{x) = $(a;) 



(38a) 



and for b(x), likewise. 



b'{x) = 
b{x) 



1 



2(mc)2 
1 



c 

e , 



2(mc)^ 

Then, by re-arranging the equation for 7, i.e., 

1 



7=1 + 6=1 



Px 



TjJ±_ Aj 

V c 



(386) 



(39) 



2(mc)^ 

an expression is provided for px in terms of the characteristics eg and Wj . 

The distribution function / = /(cq, -cuy, ro^) is an arbitrary function of its 
arguments. Now if /(eo, Wy.Wz) were to be taken as a function of eo only and not of 
voj separately then, because pj is then a function of (nj^ — eA^/c)^ as is 7, it follows 
that transverse currents 

vj^--AA (40) 
c / 

are identically zero, because the integrand is odd. Hence it is necessary and sufficient 
that f(eQ,Wy,Wz) be a function of its arguments eo and Wy and/or Wz in order to 
have spontaneous symmetry breaking and so a transverse current. 

Clearly, just as for the non-relativistic situation, different choices made for 
/(eo, OTy, -ctJz) determine the allowable soliton spatial structures. Even with the 



oc y ^ /(eo) 
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same choice of functional behavior in distribution functions for the non-relativistic 
and weakly relativistic situations one has different soliton structures. Consider, for 
example, the case of infinitely massive ions and mobile electrons with the electron 
distribution function being taken as 



fe{eo,rUy,zUz) =/oexp( -— ) exp 



fo = const. 



(41) 



Now in the weakly relativistic limit one can write the characteristic constant, eo, as 

1 



Co 



-b{x) 



l + e^(") / dx' a'{x')e~'''-'-'"> + 



Px 



c 



2(mc)^ 

thereby expressing px in terms of the characteristics, i. e, eo as defined in Eq. ((37l 
Now consider the transverse current integral 



, (42) 



_ e f d^p f e \ 

m J ■y \ c / 



(43 a) 
(43&) 



Using 1/7 1/(1 + e) f=a 1 — e w exp(— e), one can then write, in the weakly relativistic 
limit. 



J± 



e 

m 
X exp 



dpx dwy dwz 



exp 



/o 



1 



eo 

2(toc)^ 



c 



(44) 



With eo being expressed in terms of p^ and {■^^x — eA±/c)'^ through Eq. (j42|) the 
integral over px can be performed immediately yielding 

b{x) \ p p-b{x) / rx N 

l + e''(^) / da;' a'(x')e-''(^') 



J± 



2tt 1 



exp 



fo 



dwy dzUz I 



Aj_ ] exp 



X exp 



2{mcf 



A±_ 

c 



,-bix) 



(45) 



The double integral over -ojy and Wz can be done in closed form so that one can write 



w, 







X exp 



(l + « 

.2 a 



exp 



-b{x) 



& l + /3w^ 



(4 + A 



da;' a'(a;')e-''(^ ^ 



/o 



(46) 



Even in this seemingly simple extension of the non-relativistic results to the 
weakly relativistic situation one is challenged by an exceedingly non-linear set of 
equations for the field components. How many solutions the equations admit, how 
the elctrostatic and electromagnetic components are coupled, and how the solution(s) 
structure depends on the various parameters remain analytically intractable but are 
likely best addressed by using numerical procedures. 

The point being illustrated here is that, despite the characteristics being available 
in closed analytic form, the non-linear complexities of the current distribution as 
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so that 



functions of the electrostatic and electromagnetic potentials are less than inviting. 
And for each choice of distribution function similarly complex results obtain. It would 
seem that only numerical procedures can help. 

4.2. Fully relativistic behavior 

Because the characteristic Eq. ((33|) is not solvable analytically in general, two limiting 
cases will be considered in what follows. 

4-2.1. The case \a'\ <C \b' \ /j Here write the characteristic equation as 

7^ +6' = -a'7. (47) 

Then 

- 7^ + 6 = - / d^ a'7 + A = -07 + / d^ a-/' + X. (49) 

In the integral on the right-hand side set, on evaluating 7^/2 + 6 ss A, 

7' = , ^' (50) 

V2(A-&) 

to lowest order. Then 

^ 7^ + 5 + a7 + / dC = = A (51) 

2 ' J ^2 (A - b) 

and / is a function solely of A, ujy^ and vjz. 

4.2.2. The case |a'| > |6'| /7 Here write 

l^+a'.-^ (52) 
ox 7 

so that 

±{j + a) = -- (53) 
ax 7 

with, then, 

r b' 

7 + a = -y dC- + A. (54) 

In the integration on the right-hand side set 7 = A — a, (i. e., solve for 7 by neglecting 
the integral itself) to obtain 

r b' 

7 + a + / de = A (55) 

J A - a 

to lowest order, which defines the characteristic to order |6'/(a'7)| as required. Then 
/ is a function solely of A, njy, and — in contrast to Sec. 14.2.11 — of a. Because \a'\ is 
large, a direct dependence on a is the most useful choice in this casejjj] 



II Because / is an arbitrary (but positive) function of its arguments one can choose many different 
functional forms to illuminate particular points, which is what is done throughout the paper. 
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Note that the characteristic Eq. ([55]) is symmetric in {-c^7± — eA±/c)'^ so that if the 
distribution function is not exphcitly dependent on Wy,Wz but only on the constant 
arising from the characteristic Eq. (1331) . then A± = by symmetry arguments. Thus, 
just as for the non-relativistic limit: only when the distribution function depends 
explicitly on -ojy and/or Wz, in addition to the constant arising from the characteristic 
equation is there the possibility of a self-consistent electromagnetic soliton because of 
the asymmetry of the current integrals with respect to t7?j_ — eA±/c. However, the 
functional form of such solitons (and also the modifications brought about by coupling 
of the electrostatic and electromagnetic fields) is considerably different than their non- 
relativistic counterparts due to the relativistic limitation that particle speeds must be 
less than c. 

To illustrate this basic point consider again an electron-positron plasma in which 
the two distribution functions are identical. To make the comparison as close as 
possible between the relativistic and non-relativistic situations consider the electron 
and positron distribution functions [as functions of -oDy^tUz, and A, the characteristic 
from Eq. ((54)) ] to be given by 



/ = /o exp 



(56) 



with /o, Ao, and njQ constants. 

Then consider the current integral 



Jy 



(my - - Ay) f. (57) 
7 V c / 

While the general relativistic equation is solvable analytically only in the non- 
relativistic and weakly relativistic situations, it has the property that it depends only 
on the combination {u7± — eA±/c)^ . Thus the characteristic constant for the equation 
(say. A) is also a function solely of {'c^7± — eA±/c)'^. 

Hence, for distribution functions that are functions in the form /(A, njj,, cuj.) 
it follows that if / is chosen to be a function solely of A and not of ■cjy,Wz then 
all transverse current components are precisely zero. Under such conditions there 
are no electromagnetic soliton solutions. This aspect has already been seen in the 
non-relativistic (Sec. [3]) and weakly relativistic (Sec. 14. ip cases, and is now of general 
validity. Thus / must be a function of Wy and/or Wz as well as A in order to obtain an 
electromagnetic soliton. Despite the fact one cannot solve the characteristic equation 
in closed form for a fully relativistic plasma as shown in [Appendix B[ one can obtain 
accurate approximate solutions when the electrostatic field is either small or large 
compared to the Lorentz force per unit charge. In both cases it is then possible to 
express 7 in terms of the characteristic constant, as also detailed in [Appendix B[ 

The existence and structure of any solitons (electromagnetic and/or electrostatic) 
then depends on the choices made for the particle distribution functions, as also 
exhibited in detail for the non-relativistic and weakly relativistic solutions. 



5. Summary and Discussion 



While the linear Weibel instability has been thoroughly investigated, regrettably 
the same cannot be said of the non-linear behavior including the coupling of 
electromagnetic and electrostatic effects. The exploration of the non-linear aspects 
given here has uncovered a variety of effects that are germane to future investigation for 
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both non-relativistic and relativistic plasmas. While in both the non-relativistic and 
the relativistic case, the classical constants of motion are total energy and generalized 
momentum, the problem is to obtain a closed form expression for 7 (or, in the notation 
used here, Px) in terms of these constants. Because of the coupled effects of the 
magnetic fields and the electrostatic fields such appears not to be an easily tractible 
problem as shown in text. 

It has been shown previously that Weibel isolated modes (which, subsequently, 
can develop soliton modes) are retained in analytical calculations even if one allows 
for "classic" extended unstable wavenumber ranges. Because such structures develop 
only in asymmetric plasmas, they serve as an indicator for the asymmetry of the 
particle distribution function. Therefore, because precisely symmetric plasmas are 
difficult to achieve in natu r e, i solated Weibel modes will be ubiquitous, as shown 



are allowed to have a "weak" propagating component, the isolated Weibel modes are 
still generated. Hence, soliton structures should always be taken into consideration 
when investigating: (i) instabilities in (relativistic) plasmas in general; (ii) non-linear 
behavior of the resulting unstable modes; (iii) particle radiation patterns due to 
scattering in such modes. 

Perhaps the most significant theme is that the occurrence of a non-linear Weibel- 
like soliton requires that the distribution functions be dependent on all three of the 
characteristic constants. Without such a dependence (and in particular with no 
explicit dependence on characteristic constants perpendicular to the spatial variation 
direction of the soliton) then there is no electromagnetic current and so no soliton. This 
point was demonstrated for non-relativistic, weakly relativistic, and fully relativistic 
plasma situations. 

Even then, the functional behavior of the distribution functions on the three 
characteristic variables was shown, by explicit examples, to play a fundamental role in 
determining the structure of the non-linear equations for the coupled electromagnetic 
and electrostatic fields. Cases were given where no soliton was possible, where solitons 
existed only for decoupled electromagnetic field with no electrostatic component, and 
where changes in the distribution functions altered the non-linear field equations 
so markedly that each situation had to be considered anew. The characteristic 
constants could be written down in closed form for the non-relativistic and weakly 
relativistic situations, and the constants could be approximated in the fully relativistic 
plasma situation for weak electrostatic (electromagnetic) effects on an electromagnetic 
(electrostatic) field. 

Nevertheless the complexity of the resulting non-linear field equations is daunting. 
Except for simple situations it has so far not proven possible to solve such non- 
linear equations in either particular cases or the general case for chosen distribution 
functions. One suspects that only numerical procedures will allow deeper insight into 
the classes of functional behavior for distribution functions that allow solitons, for the 
spatial structure of such solitons, and for the relative contributions of the electrostatic 
and electromagnetic fields to any such solitons. Future work should attempt to 
investigate the modifications of the radiation pattern due to particle scattering in 
such soliton structures. In doing so, the question can be explored if and to what 
extent the radiation spectrum in relativistic outflows deviates from pure synchrotron 
radiation. 



recently 




Furthermore, even if the unstable wave modes 
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Appendix A. Non-Relativistic Integral Transformations 

Transformation of the non-relativistic charge and current density integrals proceeds 
as follows. In terms of the variables 

w =pl + 2m^ea^+{zUy-—Ay\^+(w^-—A^\ (A.l) 



^(y,z) = P(y,z) H A(y,z) (A. 2) 



c 



one has 



and 



Iq = j <3?p f{w,Wy,m,,) (A.3) 



c 



& p [-uj^- — AA f{w,Wy,Wz), (A.4) 



where -c^r^ = (0, Wy, vuz) and A± = (0, Ay, A^). The volume element d^p is written in 
terms of dw dwy dwz using a conventional Jacobian procedure, yielding 

d?p = „ , ^ A.5 

with Wi, = 2eaTOQ,$ + — CaA^j c)' . 

Then, by taking into account that according to the definition Eq. (jA.l|) . w is 
symmetric in -^Px and that hence one has two integration regimes Wi, ^ w < oo for 
the integration range —oo<px<oo, one obtains 

/•oo /"OO /"OO 

Iq = dw druy / dvjz {w - Wi,y^^'^ f{w,TUy,TUz), (A.6) 

Jw^ J —oo J —oo 

By defining ^ — w — w^,, this integral can be expressed as 

/•oo poo poo Q 

= -'^ / dtuj^ / dwz y/i-^ fiw^ -\- ^,Wy,n7z) . (A.7) 

Jq J-oc J -oo Cl^ 

Likewise, the current integral from Eq. (jA.4|) takes on the form 

/•oo (*00 /'OO 

Ij^ - 2 dU dwy / dwz 

■^0 J —oo 'J —oo 

X (^■:^^-^AA_'j-^f{wi,+^,zUy,mz). (A.8) 
Eqs. (|A.7|) and (IA.8I) are used in text. 
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As for the nonrelativistic situation, transformations in the relativistic situation rely 
upon the connection between the momentum variables and the canonical variables 
X^Tziy, and zuz- A simple Jacobian transformation shows that 



where one has 



with 



dx 



1 (mc) 



2 Ip.I 

-a'(x) 



dA dvoy dzuz , 



■b'{x) 



+ b{x) + / dx' a [x)^{x') = A. 



From the derivative of Eq. (jB.Sp with respect to A one obtains 



1 



1 

2 ~d\ 
or, from Eq. (|B.2p 

d_ /d-f 

Then, generally, one can write 

dX 

Note also that 

aA^^'P 



dx' 



I a'{x') dj{x 



l\2 



-fix') d\ 



exp 



a' (a;) d^'^ 
7 d\ 



dx'^^ 
j{x') 

,6V) 
7(x')2 



> 0. 



dx' 



> 0. 



(B.l) 
(B.2) 
(B.3) 

(B.4) 
(B.5) 
(B.6) 
(B.7) 



The representations from Eqs. (jB.6p and (jB.7p are particularly useful when discussing 
the limiting cases ^ \a'j\ and ^ I^'tI- 
Because 7^1 always it follows that 



A ^ i + b'{x) + I dx' a'{x')-/{x') = Aix) 



for all X so that A(x) must have a maximum value. 

Now Eq. (IB. 21) is not solvable analytically in closed form for arbitrary a{x) and 
b{x). Thus it is appropriate to investigate the limiting cases of |6'| S> \a'j\ and 
l&'l ^ \a'j\. Consider each in turn. 

Appendix B.l. Case 1. > |a'7| 

Here terms to leading order in a' will be retained because they represent the small 
order electrostatic corrections to an otherwise pure electromagnetic situation. One 
can also write 

, b'{x') 



7 + a{x) + 



dx' 



7(x') 



(B.9) 



which, with 7^1 everywhere, implies the constant v satisfies v ^ Vm, where Vm is 
the maximum value of the left-hand side of Eq. (|B.9p . 
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A connection between A and j/ is obtained simply by inserting 7 from Eq. (jB.9p 
into Eq. ((R8)) so that of Eq. (lR3t 

2 



ly — a{x) — / dx 



Consider then the volume element 



+ h{x) + / dx' a'{x')-t{x') = A. 
^0 



(B.IO) 



1 {mcY 



2 \Px\ 
{mcY 

\Px\ 



dX 



exp 



dA dzDy dwz 



\x' 



7(0;') 



dA dvoy dwz ■ 



(B.ll) 



Write 



\Px\ 

mc 



7^-1 



c 



1/2 



2 ( A - b{x) - I dx' a'{x')-i{x') 



1 



e x2 

A 

c 



1/2 



^^„^ . (B.12) 

Interest centers on the expansion of the volume element to lowest order in a{x). Then 



and 



so that 



dx' , , _ 
dx' a'{x')^{x' 



dx' a'{x') Uj2 {\~h(x')\ 



dx' c^{x')^2[\~h{x')\ 

( I a'(x') 

d^p ~ (mc) dA dtUj^ dnj^ 1~ / da;'- 
y Jo 



(B.13) 
(B.14) 



V2 [A - b{x')] 



X <^ 2 [A - b{x)] - 1 



{mcy 



e ^ 2 



-1/2 



X {1+ dx' a'{x')^/2[X-b{x')] 



2 [A - bix)] - 1 - ( 









-'1 






y] 






c 







(B.15) 



Collecting terms to leading order in a yields 
d'^p ~ (mc) dA dwy dvo^ < 2 [A — b{x)] — 1 



e s2 
•cir^ A 

c 



-1/2 



1_ /"dx' , + rdx' a'ix')W2 [X-b{x')] 

Jo y/2 [A - b{xO] Jo 



2 [A - 6(a;)] - 1 









-1 












c 







(B.16) 
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In this situation, representing an electrostatic field dominating over electromagnetic 
contributions one retains terms to leading order in b. Then one uses again 



97 



exp 



7(2;')= 



and 



so that 



Hence 



/ " w b'{x') 



exp 



7(0;')' 



3^ xjncY 97 



A-'p 



\Px 



■7 



dv 



(B.17) 
(B.18) 
(B.19) 



To leading terms in b one can replace 7(2;')^ in Eq. (jB.18[) by [ly — a(x')]^; equally 
one can write 



V — aix) — dx — 

Jo V- a{x' 



Then noting that 



\Px\ = mc 



7^-1 



one has 
d^p — {mc) 



V — a{x) — I dx 



, b'{x') 



g \2' 
T^J_ A_L 



-,1/2 



(B.20) 



(B.21) 



iy-a{x') 



Jo [v — 



V — a{x) — / da; 



, b\x') 



v~a{x') 
-1/2 



1 



{mcy 



e n2 

-1X7^ A 

c 



di/ dwy dvjz- 



(B.22) 



Expanding the right-hand side of Eq. (jB.22l) to leading terms in 6 yields 

.21-1/2 

d?p — {mc) dv dwy dwz 



2 1 

{v — a{x)) — 1 



(rnc)^ 



e ^2 

-1X7^ A 

c 



X [ly- a{x)] <l-[iy- a{x)] / dx' 

Jo v-a{x') 

dx' a{x)] dx' ^ ' 



[ly ~ a{x')Y 
[v — a(a;))^ — 1 



V- a{x') 
-1 ~i 



e x2 
T^±_ A±_ 



(B.23) 



In addition to the volume element d^p transformed to the canonical variables, 
in the integrals for transverse current contributions one also requires the elements 
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d^p/j. These elements follow directly from the above discussion because in the case 
\b'\ < |a77| one has 

b'{x') 



7 ~ — a{x) 



dx' 



V — a{x') ' 



(B.24) 



while in the case ^ W /j\ one has 

7 = 





X- 


b{x) - 1 






Jo 




X- 


b{x) - / 






Jo 



1/2 



Then 



1 



1 



1 



1 



dx' a'{x')yj2[X~b{x')] 
b'{x') 



1/2 



dx' 



v — a{x) Jq V — a(x') 



7 V — a{x) 

in the case \b'\ ^ |a'/7|, while when |6'| 3> W /l\ one has 
1 1 



|l + i [A-5(x)]"^^ dx' a'{x')^2[X~b{x')] 



(B.25) 
(B.26) 

- , , ^.".j, (B.27) 

7 ^2 [A - b{x)] ' " ^ - / ^ - I 

so that one has immediate expressions by combining Eqs. (jB.26|) with (jB.23p (in the 
case \b'\ > |a'7|). 
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